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We study the dynamics of multi-junction switching (MJS): several intrinsic Josephson junctions (IJJs) in 
an array switch to the finite voltage state simultaneously. The number of multi-switching junctions (N) 
was successfully tuned by changing the load resistance serially connected to an E^Sri^Lao^CuOg+a I J J 
array. The independence of the escape rates of N in the macroscopic quantum tunneling regime indicates 
that MJS is a successive switching process rather than a collective process. The origin of MJS is 
explained by the gradient of a load curve and the relative magnitudes of the switching currents of 
quasiparticle branches in the current-voltage plane. 



Much work is focused on the switching dynamics of current-biased Josephson junctions 
(JJs) made of cuprate superconductors. An intrinsic Josephson junction (UJ) that utilizes the 
c-axis transport of cuprates 1 ) is the most promising candidate for realizing a high-T c qubit. 2,3 ^ 
It is well known that the advantages of the IJJ are its high critical temperature and the 
excellent quality of the insulator. On the other hand, the appearance of distinctive switching 
dynamics is expected owing to the coupling effect between adjacent IJJs in an array, because 
IJJs are separated by only a few atomic layers. Therefore clarifying the switching dynamics 
in IJJs under the influence of the coupling effect is an important issue for realizing a high-T c 
qubit. Reflecting the configuration, the current-voltage (I-V) curve of an IJJ array exhibits 
a multiple-branch structure [see Fig. 2(a)]. Here we denote the critical current and switching 
current of the jth branch as I c j and I sw j, respectively (j=0,l,2,..). 

Two types of distinctive dynamics have been reported for Bi2Sr2CaCu20g +< 5 (Bi2212) IJJs 
thus far. One is a successive switching process, in which each single junction switching (SJS) 
occurs at different times when the bias current is gradually increased. This process has been 
observed in IJJs of an array having nonuniform critical currents (I c s). The properties of the 
first switching (0th branch to 1st branch) are mostly consistent with those of conventional 
J J switching: switching due to the thermal activation above the crossover temperature (T*) 
and due to the macroscopic quantum tunneling (MQT) below However, for higher- 

order switching, that is, the switching from the jth branch (j >1), the switching properties 



1/9 



Submitted to Applied Physics Express 



start pulse 




Fig. 1. (a) Schematic of the electric circuit used in the measurements, (b) Schematic of the time 
sequence illustrating the bias voltage (V&) and the junction voltage (Vjj). The start pulse is generated 
when Vb exceeds mV. The stop pulse is generated when the junction voltage exceeds Vthi- 



are significantly modified by the effect of the IJJs in the finite voltage state. 4 ~ 6 ^ A large 
enhancement of the fluctuation, detected as an increase in the escape temperature, seriously 
suppresses the switching current {I S wj-> j >1 ) as compared to the conventional J J switching 
cases. The other distinctive dynamics is a collective switching process, referred to as uniform 
junction switching, in which all IJJs in an array switch simultaneously without staying in 
any intermediate stable states. This process has been observed in IJJs with uniform values 
of I C V The escape rate of the collective switching shows TV 2 times as large as that of SJS in 
the quantum regime. Here, N represents the number of junctions switching simultaneously. 
The significant enhancement of the escape rate indicates the presence of peculiar tunneling 
processes. 8 ~ 10 ) Also the mechanism why the intermediate stable states are skipped in this 
uniform junction switching is still unclear. 

In this paper, we report on a novel type of switching process, in which the simultaneous 
switching of not all but only several of the IJJs in an array occurs. We refer to this process 
as multi-junction switching (MJS). We found that the number of simultaneously switching 
junctions N in MJS can be tuned by changing the load resistance (Rl) serially connected to 
the IJJs in the measurement circuit. The switching current distribution and the escape rate 
of MJS are systematically measured as functions of N to reveal the differences among various 
switching processes. 

We used Bi2Sri.6Lao.4Cu06+5 (Bi2201) IJJs for the present study because both the 
nonequilibrium effect and the coupling strength due to the electromagnetic effect are en- 
hanced for single-layered superconductors. The Bi2201 IJJs were fabricated by a focused ion 
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beam (FIB) etching process from a single crystal grown by the floating zone method. Details 
of the fabrication process of Bi2201 IJJs are described in ref. 11. A schematic of the electric 
circuit is shown in Fig. 1, and details of the measurement system are described in ref. 12. The 
bias current was introduced from the voltage source to the junction via a load resistance of Rl- 
A trigger circuit generated a start pulse when the bias voltage (V5), which was ramped up at 
a constant rate, exceeded Vthi- The amplified junction voltage (Vjj) was forwarded to a com- 
parator (CP), which generated a stop pulse when the output signal exceeded Vth2 as shown in 
Fig. 1(b). The time interval between the start and stop pulses was measured by a time-interval 
counter. The time distribution can be converted to a switching current distribution 4 ' 13 ) shown 
as Fig. 2(b). Typical sweep rate was O.OllA/sec and typical number of samplings to obtain 
the distribution was 10000. For the measurements of higher-order switching processes, Vth2 
was tuned to an appropriate value. 

Figure 2(a) shows typical I-V curves of a Bi2201 IJJs for Rl=500 Q and 1.2 k$7 obtained 
at T=40 mK. It is important to note that two curves are obtained for the same sample. The 
I-V curves appear to correspond to N=l for Rl=500 Q and N=2 for 1.2 kS7. To remove the 
ambiguity of N, we measured the average time (T ave ) and standard deviation (<tt) of the time 
interval distribution as functions of Vfh2- T ave and ot obtained for different Rl values are 
shown in Figs. 2(c) and 2(d). According to our previous study on SJS, 4 ) a large enhancement 
of ot is expected for higher-order switching processes because of the enhanced fluctuation. 
Actually, when Rl=500 fi, T ave has a time lag (AT) and gt is enhanced when Vth2 exceeds 4 
mV as shown in Figs. 2(c) and 2(d). This result indicates that the first switching (0th branch 
to 1st branch) and the second switching (1st branch to 2nd branch) occur separately and 
that there is an intermediate stable state on the first branch [see Fig. 3(c), LCi]. Note that 
AT of Fig. 2(c) also gives information of the actual time staying in the quasiparticle branch. 
In contrast, when Rl=1.2 kS7, neither T ave nor ut change, even for Vth2 >4 mV, within an 
accuracy of lxlO -7 s. The simultaneous switching without staying intermediate states is a 
manifestation of MJS corresponding to switching from 0th branch to 2nd branch [see Fig. 3(c), 
LC2]. It is interesting that the number of simultaneously switching junctions can be tuned 
simply by changing the magnitude of Rl- 

Figure 3(a) shows a schematic diagram of N as a function of Rl and the temperature deter- 
mined by the above-mentioned procedure. Symbols in the figure correspond to experimentally 
detected points. It is clear that N increases with increasing Rl for a fixed temperature. Taking 
account of this feature, the origin of the MJS process can be phenomenologically explained as 
follows. The most important point is the relative magnitudes of the I c s and I sw s, as schemat- 
ically shown in Fig. 3(b). I c q is always smaller than I c \ and I C 2 (I c o < Ici,Ic2) because each 
junction switches to the finite voltage state in order of increasing I c when the bias voltage is 
ramped up. In contrast, since the escape temperature for higher-order switching processes is 
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markedly increased, I sw \ and I SW 2 are considerably reduced from I c \ and I C 2, respectively. 4 ) 
This enhanced escape temperature can induce the reversal of the orders of the magnitudes, 
giving I sw q > I sw i, I S w2- The states on the first and second branches are no longer stable for 
I c i > I > I sw i and I C 2 > I > I S w2, respectively. The MJS process is promoted by these " 
missing quasiparticle branches (MQBs)". Figure 3(c) schematically illustrates the trajectories 
of the SJS and the MJS processes while ramping up the bias voltage. When the bias voltage 
is gradually increased, the junction state jumps from A on the Oth branch to Bi or B2 on 
the 1st branch along the load curve whose gradient is inversely proportional to Rl- The line 
connecting I sw o and I sw i, labeled as LCo (693.2 £1 in the present IJJ) in Fig. 3(c), corresponds 
to the critical load curve that separates SJS and MJS. When Rl < Rlo, the junction state 
jumps from A to Bi on the stable quasiparticle branch; therefore MJS does not occur. In 
contrast, when Rl is considerably larger than Rlo, the state jumps from A to B2 that locates 
inside the MQB, then the next jump immediately occurs without staying on the branch. This 
immediate switching is the origin of MJS. Actually, the critical load curve obtained in this 
scheme is comparable to experimentally obtained value. Meanwhile, when Rl is comparable 
to Rlo, more complex dynamics is expected, which will be described in a separate paper. 

The above analysis suggests that the MJS process originates from the stability of higher- 
order branches rather than the onset of the collective switching process. To confirm this 
scheme, we measured the temperature dependence of the standard deviation of the switching 
current distribution P(I) for Rl=500 $7, 1.2 kSl and 5.2 kO as shown in Fig. 4. Vth2 was 
fixed at 1 mV. For all cases, the switching current distribution forms a single curve that has 
T 2 / 3 dependence above T* ~0.23 K dominated by the thermal activation and is saturated 
below T* dominated by the MQT. Theoretically calculated T* of 0.29 K based on eq. 19 of 
ref. 13 using the plasma frequency of 88GHz, the capacitance of 212 fF, and quality factor 
of 22 obtained by microwave irradiation is consistent with the experimental value. These 
results confirm that the switching features of MJS mostly agree with those of conventional 
JJs. The A^ 2 enhancement of the escape rate expected for collective switching cannot be 
detected at present. Therefore, we conclude that each switching in the MJS process takes place 
successively rather than collectively. However, MJS is distinct from conventional successive 
switching in the point that the enhancements of the fluctuation in the higher order switching 
processes are completely absent as previously mentioned. 

An important question is why the collective switching process reported in ref. 7 cannot be 
observed in the present measurements. We speculate that this inconsistency stems from the 
difference in the uniformity of the values of I c in IJJs in an array. The values of I c in an IJJ 
array fabricated by the FIB process tend to be more widely distributed than those fabricated 
by Ar ion milling. 7 ) Actually, the estimated values of I c q and I c \ based on the Kramers escape 
rate and the MQT models are 21.2 and 21.5 uA, respectively. The present experimental 
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results verify that these values are sufficiently close for the emergence of MJS, but probably 
not sufficiently close for the emergence of the collective switching process. Therefore, it will 
be interesting to clarify how much uniformity is required to observe the collective process. 
This issue will be confirmed in future works. 

Acknowledgments This work was financially supported by a Grant-in- Aid for Young Scientists 
(b)(No. 21710100) from Japan Society for the Promotion of Science, Japan, and by Mitsubishi foun- 
dation. 



5/9 



Submitted to Applied Physics Express 



(a) 



(b) 




time [ms] 

1.9 1.91 1.92 1.93 



(c) 



10 15 
V [ mV] 



20 




experiment 
■ calculation 



8 



„ 6 

> 

EM 

> 2- 




, AT . 



X X 500Q 
K O 1.2kQ 



(d) 



20.4 20.5 20.6 20.7 20.S 20.9 

I [MA] 

8 



LS 



2.2 



„ 6 

> 

> 2 




X 
X 



X 500Q 
O 1 .2kQ 







T ave [ms] 



20 40 



60 



Fig. 2. (a) Typical i-V curves of the Bi2201 IJJ in the cases of Rl— 500 ft and 1.2 kSl measured at 
T=40 mK. (b) Typical switching current (or time interval) distribution obtained at T=40 mK. Points 
represent experimental data. Solid lines represent the theoretical fitting at MQT region, (c) (d) Rl 
dependences of the average time (T ave ) and the standard deviation of time (or) as functions of Vthi 
obtained from the time interval distribution at T—40 mK. 
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Fig. 3. (a) Schematic diagram of ./V as a function of Rl and temperature. Symbols in the figure corre- 
spond to experimentally detected points, (b) Schematic of the I-V curve and the relative magnitude of 
the critical currents and the switching currents. Missing quasiparticle branches (MQBs) corresponding 
to the unstable states (/ > I sw i,I sw ^) are shown as dotted lines. (C) Schematic illustration of the SJS 
and the MJS processes in the I-V plane. The trajectories of SJS and MJS while ramping up the bias 
voltage are LCi and LC2, respectively. As the bias voltage is gradually increased, the state jumps from 
A on the Oth branch to Bi or B2 on the first branch along the load curve whose gradient is inversely 
proportional to Rl- For the case of LCi, SJS occurs because Bi is stable, while for the case of LC2, 
MJS occurs because B2 is on the MQB. 
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Fig. 4. Temperature dependence of 07 denned by the standard deviation of the switching probability 
distribution P(I) for Rl = 500 f2 (crosses), 1.2 kfi (circles), and 5.2 kfi (triangles). T 2 / 3 dependence 
in the thermal activation regime and saturation in the MQT regime (below 0.23 K) are detected. The 
inset shows the escape rates for N=l, 2 and 4 at 40 mK. 



8/9 



Submitted to Applied Physics Express 



References 

R. Kleiner, F. Steinmeyer, G. Kunnek, and P. MiMer: Phys. Rev. Lett. 68 (1992) 2394. 

K. Inomata, S. Sato, K. Nakajima, A. Tanaka, Y. Takano, H. B. Wang, M. Nagao, H. 
Hatano, and S. Kawabata: Phys. Rev. Lett. 95 (2005) 107005. 

T. Bauch, F. Lombardi, F. Tafuri, A. Barone, G. Rotoli, P. Delsing, and T. Claeson: Phys. 
Rev. Lett. 94 (2005) 087003. 

H. Kashiwaya, T. Matsumoto, H. Shibata, S. Kashiwaya, H. Eisaki, Y. Yoshida, S. Kawa- 
bata, and Y. Tanaka: J. Phys. Soc. Jpn. 77 (2008) 104708. 

K. Ota, K. Hamada, R. Takemura, M. Ohmaki, T. Machi, K. Tanabe, M. Suzuki, A. 
Maeda, and H. Kitano: Phys. Rev. B. 79 (2009) 134505. 

P. A. Warburton, S. Saleem, J. C. Fenton, M. Korsah, and C. R. M. Grovenor: Phys. Rev. 
Lett. 103 (2009) 217002. 

X. Y. Jin, J. Lisenfeld, Y. Koval, A. Lukashenko, A. V. Ustinov, and P. Miiller: Phys. Rev. 
Lett. 96 (2006) 177003. 

M. V. Fistul: Phys. Rev. B. 75 (2007) 014502. 

S. Savel'ev, A. O. Sboychakov, A. L. Rakhmanov and F. Nori: Phys. Rev. B. 77 (2009) 
014509. 

T. Koyama and M. Machida: J. Phys.: Conf. Ser. 150 (2009) 052127. 

T. Matsumoto, H. Kashiwaya, H. Shibata, H. Eisaki, Y. Yoshida, and S. Kashiwaya: Appl. 
Phys. Express. 1 (2008) 101701. 

H. Kashiwaya, T. Matsumoto, H. Shibata, K. Tani, and S. Kashiwaya: IEICE Trans. 
Electron E90-C (2007) 605. 

A. Wallraff, A. Lukashenko, C. Coqui, A. Kemp, T. Duty, and A. V. Ustinov: Rev. Sci. 
Instrum. 74 (2003) 3740. 



9/9 



